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Inflation typically predicts a quasi scale-invariant spectrum of gravitational waves. In models of 
slow-roll inflation, the amplitude of such a background is too small to allow direct detection without 
a dedicated space-based experiment such as the proposed BBO or DECIGO. In this paper we note 
that particle production during inflation can generate a feature in the spectrum of primordial grav- 
itational waves. We discuss the possibility that such a feature might be detected by ground-based 
laser interferometers such as Advanced LIGO and Advanced Virgo, which will become operational 
in the next few years. We also discuss the prospects of detection by a space interferometer like LISA. 
We first study gravitational waves induced by nonperturbative, explosive particle production during 
inflation, finding that explosive production of scalar or vector quanta does not generate a significant 
bump in the primordial tensor spectrum. On the other hand, chiral gravitational waves produced by 
electromagnetic fields amplified by an axion-like inflaton could be detectable by Advanced LIGO. 
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1. INTRODUCTION. 

Once we consider perturbations on the top of 
a homogeneous and isotropic Friedmann-Robsertson- 
Walker Universe, inflation generically predicts quasi- 
scale-invariant spectra of scalar and tensor perturba- 
tions. 

The scalar perturbations have been detected, and all 
their properties appear to wonderfully agree with the 
predictions of the simplest models of inflation. The ten- 
sor modes, however, have not yet been detected, and we 
put our best hopes to find them in the study of their ef- 
fect on the polarization of the Cosmic Microwave Back- 
ground photons. A direct detection of the tensor modes 
from inflation, on the other hand, is not expected before 
dedicated space-based interferometers such as the pro- 
posed BBO or DECIGO [1], are launched in the next 
few decades. Indeed, gravitational interactions are so 
weak that we have not yet detected gravitational waves 
of any origin. 

Experiments searching for gravitational waves of as- 
trophysical origin, such as LIGO [2], GEO600 [3], 
Virgo [4], and TAMA300 [5], have been taking data 
for several years. LIGO and Virgo will see their sen- 
sitivity improved by a factor of <~ 10 in the next few 
years and might detect the first gravitational wave as 
early as 2015. They will be sensitive to a stochastic 
background of gravitational waves whose logarithmic 
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contribution to the critical density, Qgw h 2 , is of the 
order of 10~ 9 at a frequency of ~ 100 Hz. LGCT [6] 
will have a comparable sensitivity at similar frequencies. 
A space-based experiment like LISA would be able to 
reach flaw h 2 ~ 10~ n at 10~ 3 Hz. The proposed Ein- 
stein Telescope [7] would have similar sensitivity while 
working at LIGO frequencies. Non-detection of tensor 
modes at CMB scales strongly constrains a scale in- 
variant background of tensors of inflationary origin to 
Slcff h 2 < 10~ 14 . Since the spectrum of primordial ten- 
sor modes is generically flat or slightly red, none of these 
detectors is usually expected to be able to detect tensor 
modes produced during inflation [8]. 

In this paper we note that production of particles dur- 
ing inflation generates a feature in the tensor spectrum 
which could be detectable by gravitational interferom- 
eters in the (relatively) near future, without conflicting 
with CMB constraints. 

Nonperturbative production of particles during infla- 
tion, first studied in [9], is possible because the rolling 
inflaton provides a time-dependent background that 
carries the energy necessary for the production of light 
species. The simplest and most studied example of such 
a system is given by a scalar field \ whose mass depends 
on the inflaton (f>. If \ becomes effectively massless as 
the inflaton rolls down its potential, then it becomes 
energetically cheap to produce its quanta. In this case, 
particle production happens at a precise moment during 
inflation determined by the time when the total mass 
of x crosses zero. A second possibility is that the in- 
flaton <fi couples to a derivative of some field such as a 
gauge field [10]. In this case, the field can stay mass- 
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less as the inflaton rolls down its potential, and particle 
production can happen steadily during inflation. 

Particles produced through these mechanisms carry 
energy- momentum tensor T^, which perturbs the FRW 
metric into: 

= a(T) 2 (-dr 2 + (Sij + hij) dx'dx 1 ) , (1) 

where we use conformal time 1 r and a transverse, trace- 
less decomposition, ignoring perturbations which do not 
source gravitational waves. The equation of motion of 
the perturbations is: 



K 3 



where Hij 



ti^-Aha = — n„ ab T ob 



M 2 13 



(2) 



II im is the transverse, 
traceless projector, and = Sij — didj/A. For any 
given expression of T a b(x, t), one can then solve for- 
mally cq. (2) as 

hij(k, t) = ^ j dr f G k (T, r')IL/ b (k)T ab (k, r') , 

(3) 

where G^{t, r') is the retarded propagator solving the 
homogeneous transform of cq. (2). In this paper we 
will assume a de Sitter background a(r) = — (Ht)~ , 
so that the retarded propagator reads 



G k (T,T>) = 



fc 3 T> 



(l + fc 2 rr') sin fc(r-r') 



+ k (t' -t) cos fc(r-r') 6 (r - r') . (4) 

We will examine several ways of generating a nonvan- 
ishing T a b in eq. (3), some of which induce a significant 
feature in the spectrum of tensor modes. If particle 
production happens explosively at a precise time dur- 
ing inflation (as in [9]), the spectrum of gravitational 
waves will show a feature at the scale corresponding 
to the time of particle production. If particle produc- 
tion happens continuously, on the other hand, then the 
spectrum of induced tensor modes will be smoother. 
It is worth noting that, since the source of the grav- 
itational waves is quadratic in a gaussian field, hij 
is expected to have a maximal three point function 
(hhh) ~ (hh) 3 ^ 2 [11-14]. The direct detectability of 



1 Throughout the paper we will denote by a prime a derivative with 
respect to the conformal time r, and by a overdot a derivative 
with respect to the cosmological time t: ' = d/dr, ' = d/dt. 



tensor modes produced during inflation was also con- 
sidered by [15], where it was shown that, if the post- 
inflationary Universe is dominated by a fluid stiffer than 
radiation, the primordial tensors amplitude can be sig- 
nificantly enhanced. More recently, [16] has studied the 
detectability of gravitational waves produced by phase 
transitions during inflation. 

In section 2 we show that the explosive production of 
quanta of a scalar field x can only generate a tiny cor- 
rection to the background, quasi-scale invariant spec- 
trum of tensor modes. A scenario - "trapped inflation" 
- where explosive production of particles occurs sev- 
eral times per efolding of inflation has been considered 
in [17]. In this scenario, particle production slows down 
the rolling of the inflaton so that inflation can occur 
even on a (relatively) steep potential. Our analysis will 
allow, in section 2 C, the evaluation of the amplitude of 
the tensors induced by the trapping fields. 

In section 3 we study a mechanism analogous to that 
of section 2, where the scalar field \ is replaced by a vec- 
tor field A^. The amplitude of the gravitational waves 
produced in this case is equal to that found in section 
2, multiplied by a factor two that accounts for the two 
helicity modes of the vector field. 

Finally, in section 4 we discuss the case of tensor 
modes produced through a gauge field coupled to an 
axion-like inflaton (as discussed in [20, 21] and, in 
greater detail, in [22]). In this case the amplitude of 
gravitational waves can increase dramatically at smaller 
scales so that the system can obey the WMAP con- 
straint on primordial tensors at CMB scales and still 
lead to detectable tensors at scales probed by ground- 
based laser intereferometers. 



2. GRAVITATIONAL WAVES FROM SUDDEN 
PRODUCTION OF SCALARS DURING 
INFLATION. 

Several systems can lead to the production of par- 
ticles during inflation. The one that has received the 
most attention is described by the lagrangian 



C 



1 



1 



d^d^-v (<j>)- -d^x 



' 2 



(0- 



M 2 x 2 , 

(5) 

where V(<fi) is the potential supporting inflation and 
where we neglect for simplicity the self-interactions of 
the field x- If the inflaton </>(r), while slowly rolling 
down V((j)), crosses the value cf) a , then the field \ be- 
comes momentarily massless, and its quanta can be co- 
piously produced. The analysis of [23] shows that the 
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occupation number of x, shortly after </> crossed </>o, is 
given by n x (k) = exp | —7r ^ | }i where 4>q ~ d(f>/dt at 
the time </> crosses <Pq. 

In this section we compute the number of gravi- 
tons produced by these quanta of x- The spatial part 
of the stress-energy tensor for the field x is given by 
T a b = d a \dbX + $ab{- ■ ■ ), where the part proportional 
to S a b is projected away by Tlij ab - We promote the field 
x(x, t) to an operator x(x, r), which we Fourier trans- 
form, factoring one power of the scale factor a(r) for 
canonical normalization 

Plugging this decomposition into eq. (3) gives the tensor 
spectrum 

1 (' d>T^ 

(fcy(k, ^)Mk', t)) = J a pp G fc(^ T ') x 

f Ht" 

X J <^ Gfe ' (T ' r'WiK^j^x 

x J d 3 pd 3 p' Pa (k b - p b )p' c (k d - p^)x (7) 

x ( X (p, r')x(k - P, r')x(p', r")x(k' - p', r")) . 

The quantity (...) in the equation above can be re- 
duced using Wick's theorem and ignoring the discon- 
nected term proportional to <5( 3 )(k)<5( 3 )(k'). Hence we 
need only to evaluate the two-point correlators, being 
careful to avoid divergences. 

The evolution of the system occurs in three stages: 

(i) At early times, significantly before <f> reaches </>o, 
the Universe does not contain quanta of the \ held. 
No gravitational waves are produced by x during this 
epoch. 

(ii) As <p gets close to (j> , the effective mass of x, 
m x (t) = g (<j>(t) — 4>o), starts evolving nonadiabatically, 
m x > m 2 . The duration Ai nad of the epoch of nona- 

diabaticity is Ai nac i ~ {g^o)^ 1 ^ 2 , which must be much 
shorter than a Hubble time for production of quanta of 
X to be efficient. During this short epoch, the quanta 
of x, while being produced, source gravitational waves. 

(Hi) After a time of the order of Ai na( j after <j> has 
passed 4>o, the mass m x (t) evolves adiabatically again. 
Even if they are not being produced, quanta of x are 
still filling the Universe and source the tensors before 
diluting away in a few efoldings. 

In the following subsection we will study the gravita- 
tional waves produced during the epoch (in), while in 



subsection 2 B we will discuss those produced during 
the nonadiabatic period (ii). As we will see, the ten- 
sors produced during these two epochs have comparable 
amplitude. 



A. Tensor production during the adiabatic epoch. 

The main quantity we have to evaluate is 
(x(P) T ')x(li T "))- We decompose x(k) into cre- 
ation/annihilation operators as x(k, 7") = x(k, t) dk + 
X* (— k, r) d^ k , where the function x must obey the 
equation 

X"(k, r)+w(k, r) 2 x(k, r) = (8) 

with 

w(k, r) 2 = k 2 + g 2 a(r) 2 (4>(r) O ) 2 - ^ , (9) 

and dk annihilates the vacuum during period (i). We 
define the Bogolyubov coefficients a(k, r) and /3(k, r) 
via 

X(k, r) =-±= (e^-T^k, r) + e^ T -/3(k, rj) 

X'(k, r) =i ^( - e- i / T " a(k, r) + e*^ " /3(k, r)) 

(10) 

so that, in the adiabatic limit a/ <C <jJ 2 , a and /? are 
constants. This way, we can rewrite 

p—i f T ujdf p i f T cjdf 

X(k,r) = - b k(r) + ^^St_ k (r), (11) 

where we have defined the new annihilation operator 
S k (r) = a(k, r) d k + [3* (-k, r) a f _ k , (12) 

that annihilates the vacuum during period fm). 

In order to renormalize the theory, we impose that 
the operator x(p, t')x(<1, t ") within (...) is normal or- 
dered. However, we require normal ordering in terms of 
the 6k operators while using the vacuum state defined 
by the dk operators. This way we calculate the number 
of quanta of our initial particle definition existing at the 
end. Using the decomposition (12) and the commuta- 
tion relation [d q , d|J = <r 3 )(p — q), we obtain 
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<*(p, ^')x(q, r")) 



^ (3) (p + q) 

2 v /^ P (r')^p(r") 



( e i/ T T "«P /3*(- p , r')/?(-p, r") + h.c.) 
+ (e"^' "p-'-T" ^a(p, r')/T (p, r") + (/ r", h.c.)) 



(13) 



Using the expression above and Wick's theorem, eq. (7) can be written as 



(%(k, r) hij(k.', t)) 



■ k k j.^pir 



8 7r 3 If; 



k 2 



) J *(r>y 



G k (r,r') 



dr" G k (r,T") 



— ^ f dr" 

) 2 ^/a;p(r') w k _ p (r') J «(t") 2 y/u p ( T ") w k - P (-r") 



x [e^' w p/3*(| P |, r')/3(|p|, t") + e- l / T ' w pe-^ T " w pa(|p|,r')/3*(|p|,r") + (r' T",h.c.) 



(14) 



x e l -/V T '" Wk -p/3*(|k-p|, r')^(|k-p|, t") + e* f T ' ^- Pe - 4 / T " <-k- Pa (|k _ p | ; r ')/3*(|k - p|, t") + (t' r",h.c.) . 



When multiplied out, some terms in the above equa- 
tion arc rapidly oscillating and we neglect them as they 
give subdominant contribution to the integrals. At this 
point we need to evaluate the Bogolyubov coefficients 
a and (3. 



1. Evaluating ft. 

The function \ obeys cq. (8), and quanta of \ are 
produced during the short epoch of nonadiabatic evo- 
lution of uj(k, r), during which we neglect the ex- 
pansion of the Universe. In a de Sitter background 
and in the slow-roll approximation, <ft evolves approx- 
imately linearly in physical time (j>(t) = (f>o + 4>o t or 

as 4>{t) — 0o — /T^°S(n)) * n con f° rma l time, t and t 
arc defined such that cf){t = 0) = 4>(t = t ) = <j)o- The 
duration Ai na( j of the nonadiabatic epoch is determined 

by the condition m x >, m 2 ,, yielding Ai na d ~ 1/ 
It is consistent to neglect the expansion of the Universe 
if Ai na( j <C 1/H so that the validity of our analysis 
requires g ^> H 2 /\<p a \. 

Under these conditions the equation for \ during the 
nonadiabatic epoch reduces to 

X+(k 2 H 2 r 2 +g 2 ^t 2 ) X = 0, (15) 

to which we can apply the analysis of [23], obtaining, 



up to an irrelevant phase, the Bogolyubov coefficients 

oi{t > T ,k) = Vl + e-™ 2 e la * , 
(3(T>r a ,k) = e-^ 2 , _ (16) 

where we have defined k = k H To/ \J g <j>o and a K = 

Arg[r(i±|^)]+^ (l-log£). 

2. The two point function. 

After the phase of nonadiabatic evolution of m x , the 
Universe contains <~ J d 3 k\/3\ 2 /a(r) 3 quanta of \ P er 
unit volume. We are now in the position of computing 
the spectrum of gravitational waves generated by such a 
population, which is given by the integral (14), where u> 
is defined in (9) . The expression of uj can be drastically 
simplified by observing that, at the end of the nonadia- 
batic period, g (0 - <j> ) ~ g <fio Ai na d ~ (5 0o) 1/2 > H. 
As a consequence, the second term in (9) is much larger 
than the third one. Moreover, the exponential suppres- 
sion in (3 means only momenta k < (g qjo) 1 ' 2 /{H t$) con- 
tribute significantly to the integral. Noting also that, 
following the non-adiabatic period, |r| < |r |, the k 2 in 
(9) is negligible as well. 

Therefore, during the entire phase (Hi) of adia- 
batic evolution of the system, we can approximate 
uj ~ \g (0(t) - 4> )/{Ht)\, where 0(r) - <f> ^ 
— (4>o/H) log(r/r ). Using these estimates, we can 
write the integral (14) in the limit r — > (i.e., at the 
end of inflation, long after r ) as 
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(Mk) M»0> 



5( 3 )(k + k') ff 8 



4 7r 3 fc 6 g 2 <j> 2 M p 



d 3 p p 



(p-k) ; 

k 2 



(17) 



x { |/3 (p)| 2 |/3 (k - p)| 2 +Rc [a(p) a*(k - p) /3*(p) /3(k - p)] } 



T + AT nad 



sin(fcr) — fcr (fcr) 
ln(r /r) 



where Ar nad /|T | = At nad ~ H (gfo)- 1 / 2 < 1. As 
we will see in eq. (20) below, the result depends only 
logarithmically on Ai nac j so that ignorance of its exact 
value does not affect significantly the results. Next, we 
recognize that at large values of fc the two point function 
is suppressed by the factor of (sin fcr — fcr cos fcr) / A: 3 
coming from the Green's functions. Since H 2 /gcb <C 1, 
p is only suppressed after pr a > g4>/H 2 ^> 1. Therefore 
the integrand gets its main contribution from the region 
p^> k. Using these approximations, the above equation 
is simplified to 



spectrum reads 



P*(k) ~ 



2H 2 
^M 2 



1 + 4 8 x 1(T 4 ( fcr o cosfcT o -siufcrp) _ 



M 2 { H 2 



3/2 



log 2 



\kr a \ 



H 



(21) 



(fcy(k) M k ')) 



47r 3 fc 6 g 2 /j) 2 M p 



dr ■ 



r +AT na 



sin(fc-r) — fcr (fcr) 
Hto/t) 



2 e si- + e s* " j / 

2 



(18) 



After computing the integral in d 3 p, we are left with 
the simple expression 



(Mk) hiiW) 



^ 3 )(k + k') H 




d/T0 , (fc|r |) , (19) 



where 



Fe(y) 



(sin xy — xy cos xy) , 

x 5 ax 

log a; 



[(y cosy - siny) loge] . 



(20) 



We thus see that the effect of the creation of quanta 
of x is to superimpose a scale dependent contribution 
to the scale invariant spectrum of tensors generated by 
inflation. 

We next observe that </> = \/2e H Mp, where e«l 
is the slow-roll parameter. Supplying reasonable val- 
ues for H, Mp, and e allows us to find that the |fcr |- 
dependent part of the spectrum is maximized at |fcro| — 
2, where the component from particle production eval- 
uates to ~ 1.8 x 10- 4 {%t? 12 log 2 (^). Using 

the same approximation, the log 2 term gives at most a 
factor ~ 10 2 ; therefore, the correction to the standard 
result is at most of the order 10~ 2 H 1 / 2 /Mp 2 , which is 
several orders of magnitude smaller than unity. 

We thus conclude that the presence of a gas of adia- 
batically evolving scalar particles produced nonpertur- 
batively during inflation generates a tiny correction to 
the spectrum of primordial tensors. This result agrees 
with [24], where it was shown that in the Minkowsky 
limit, H — > 0, a gas of adiabatically evolving scalars 
does not generate any gravitational waves. In our case, 
since we are on an expanding background, gravitational 
waves are produced, but the effect is still small and un- 
observable. 



The two point function (19) should be added to the 
standard, quasi scale invariant contribution from infla- 
tionary gravitational waves so that the resulting power 



Let us next estimate the amount of gravitational 
waves produced during the period of nonadiabatic evo- 
lution of m x (r). 
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Tensor production during the nonadiabatic 
epoch. 



The period of nonadiabatic evolution of m x (r) lasts 
much less than a Hubble time. We will use again our 
physical time variable t = H^ 1 log(To/r). Since we are 
now looking at the period \Ht\ <C 1, we can approxi- 
mate the change of variable as r ~ To (1 — Hi). This 
implies that we can replace r' and t" by To in the in- 
tegrands of eq. (7). During this short time, the field \ 
will obey eq. (15). 



During the periods of adiabatic evolution of m x , the 
concept of a particle of x 1S weu defined, and the use of 
the Bogolyubov coefficients gives an appropriate way of 
computing the spectrum of gravitons produced by the 
gas of quanta of \. During the short epoch of nonadi- 
abatic evolution of m x , however, it is more convenient 
to switch to a different prescription. Following e.g. [25], 



we set 

(x(p, Ox(q, t")) = 5< 3 >(p + q) [ X (P, i')x*(p, t")- 

-x(p,*')xW)] (22) 

where x(p, t) is the solution, with appropriate bound- 
ary conditions, to eq. (15), and x(p, i) is the solution 
to the same equation in the adiabatic approximation: 



X (p, t) = e-^/* D ^ [( _! + i)r)] 



x(p,t) 



^/JfjT e i[%\/r) 2 +P 2 +p 2 /2 log(77/p+ v /l +') 2 IP 2 )] 



Wo) 1 / 4 



(23) 



where we have defined a dimcnsionlcss time r\ = 
(g^oY^t and a dimensionless momentum p — 
pHro/igipo) 1 an( i D a {z) is the parabolic cylinder 
function. The term in x(p, i) in eq. (22) takes care of 
the UV-divergent terms which would otherwise appear 
in the tensor spectrum. 

Working forward from eq. (22), the expression for the 
contribution of the nonadiabatic epoch to the two point 
correlator of the graviton simplifies to 



(hijCk, t) hij(k', t)) 



6^ (k + k') H 4 

2 7T 3 M% ~W 



sin k tq — k To cos k tq 



d 3 p p 



(P"k) ; 

k 2 



(Ht ) 2 



(24) 



J dt'dt" [ X (p, t') X * (P, t") - x(p, t')x* (p, *")] [X(k - P, t') X * (k - p, t") - x(k - p, t')x* (k - p, t")] 



We next observe that for \kTo\ ^> 1 the two point 
function of the graviton is suppressed by the coeffi- 
cient oc fc~ 4 in front of the integrals of eq. (24). In 
the unsuppressed regime \kT \ £ 1, the quantity k — 
kH \To\/{g4>o) 1 ^ 2 "C 1, so that we can set k = in 
the second line of eq. (24). At this point the angu- 
lar integral can be easily computed, and by observ- 
ing that both x an d X m e q- (23) are of the form 
(H ItqI) 1 / 2 (g 0o) -1 / 4 x [function of (p, if)], we can write, 



after appropriate changes of variables, 



(hij(k, t) hij(k', t)) 



16 (Jf 3 ) (k + k') H 4 



15 7T 2 fc 6 |r | 3 Mp 



H 2 



3/2 

| [sin k t — k r cos k r ] 2 x 

x Jfdp J dr( d V " F(p rf, rl'f , (25) 

where F(p, 77', rj") in the second line is a dimension- 
less function built out of the dimensionless functions 
appearing in eqs. (23), whose integral gives an 0(1) re- 
sult. 
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The result of this subsection is that the contribution 
to the two-point function of the graviton by the nonadi- 
abatic epoch (ii) has the same form, modulo a logarith- 
mic term, as the contribution from the adiabatic epoch 
discussed in section 2 A. We conclude that the overall 
effect of graviton creation by the scalars \ gives a neg- 
ligible correction to the standard spectrum of tensors 
generated by inflation. Let us next discuss what hap- 
pens if there are several events of explosive production 
of scalars. 



i.e., in an order of magnitude estimate, 



H 4 



trapped — (2 7r) 3 if A Mp \ H 2 




3/2 



(28) 



which shows that the resulting spectrum is scale invari- 
ant. 

Let us now connect this result to the parameter space 
of trapped inflation. Plugging the "slow roll" equation 
of trapped inflation [17], {g<j) ) 5/2 ^ (2 tt) 3 H A V, into 
the equation for the tensors, we obtain 



C. An application: tensor modes in trapped 
inflation. 



•p* , ~ 

' trapped — „^4 



(29) 



The analysis of sections 2 A and 2 B has shown that 
the spectrum of tensor modes induced by a single event 
of production of scalars has the form 



H A (fc|r | cos fc|r | — sin fc|r |)^ 
M% (2 7r) 3 fc 3 |r | 3 




(26) 



up to a coefficient of order one. Let us now suppose 
that there are several instances of particle production. 
In particular, we will assume that particle production 
happens so often to lead to trapped inflation [17], i.e., 
to a slowing down of the rolling of the inflaton more 
significant than that due to Hubble friction 2 . Since the 
tensor modes excited by each burst of particles sum 
incoherently, the resulting power spectrum of the tensor 
will just be the sum of the individual power spectra. If 
the bursts happen frequently enough, summation the 
various contributions to V l x is equivalent to integration 
over ^ J^-, where A = fc+i — fa is the distance in 
field space between two sites of particle production. We 
therefore get, 



V 1 

' trapped 




cos fc|ro| — sin fcl-ro I 



(2 7t) 3 /c 3 |t | 3 



(27) 



2 The mechanism leading to trapped inflation is similar to that 
responsible for warm inflation [18, 19], although in the case of 
warm inflation friction is provided by a thermal bath rather than 
by nonperturbative particle production. 



This should be compared to the standard amplitude 
of gravitational waves 'P* tandard ~ H 2 /Mp, so that 



V 1 

trapped 



/Standard - HV'/gV. In order to proceed 
further, we choose a form of the potential. For V((j>) = 

m 2 4> 2 /2 (the case Studied in [17]) trapped /^standard - 

m/gMp. This is much smaller than unity in the 
phenomcnologically allowed region of parameter space 
of [17]. Therefore, even in the case of several bursts of 
scalars, the induced spectrum of tensors is subdominant 
with respect to the standard one ~ H 2 /Mp. 



3. PRODUCTION VIA VECTORS. 

In this section we study the production of gravita- 
tional waves induced by the explosive production of 
vectors. We will only focus on the gravitational waves 
produced during the period of adiabatic evolution sub- 
sequent to the creation of the vectors. Based on the 
similarities between the results of sections 2 A and 2 
B, we expect that the contribution of the nonadiabatic 
period will at most give an order one correction to the 
results presented here. 

We consider the gauge-invariant lagrangian 



1 



F^F^ -(£>"*) (D^)*-Vm 2 ), (30) 



where = V M 



i e is the gauge-covariant deriva- 
tive. We will assume that the Higgs field \& is a function 
of time, but is otherwise homogeneous during inflation, 
Defining = ip e 10 , it is possible to show [26] that, if * 
is spatially homogeneous, it is consistent to choose the 
gauge A (t, x) = 0(t)/e, V • A = 0, with A satisfying 
the equation 

A"(k, t) + (k 2 + e 2 a 2 (r) ^ 2 (r)) A(k, r) = , (31) 
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so when ip crosses zero, production of photons occurs 
precisely in the same way as described above in the 
case of scalars. In particular, the expression of the 
Bogolyubov coefficients is calculated the same way as 
that of section 2.A.I. The portion of spatial part of the 
stress-energy tensor of A we are interested in is given by 
Ttj = A\ A'j + e m d k A l e jmn 8 m A n - to 2 (t) A t A jt where 
we have defined m(r) = ea(r) ip(r). We promote the 
photon to a quantum field A(r, x), which we decompose 
on a basis of helicity vectors e± (k) satisfying the con- 
ditions k-e±(k) = 0, e CT (k) -e CT '(k) = 5 a - a i, e±(— k) = 
— e T (k), k x e±(k) = =Fifce ± (k) and e±(k)* = e T (k): 



A=± J 



d 3 k 

(2tt) 3/2 



e A (k)i A (r, k)e tk x + h.c. , (32) 



with A\(t, k) = Ax(t, k) o(k). Eq. (3) then reads 



hx(K t) = 



Ml 



E / dT 



, G k (T, T 1 ) 



d 3 q 



a(r') 2 J (2tt) 3 / 2 
x nf(k) e l M) e£(k - q) [i'^q, r') i' a ,(k - q, r') 
+ {a a' q |k - q| - m 2 (r')) i ff (q, r') A a , (k - q, r') 



(33) 



where we have introduced the projector II A (k) = 
e A (k) e A n (k)/\/2 on the helicity- A component of the 
graviton. 

We are now in position to compute the two-point 
function of the graviton (/i A (k, r)/i A '(k', r)). To do 
so, we use the Wick theorem to express it as a function 
of two-point correlators of the gauge field. Next, we 
define a function /2(k, r', r") such that 



(i ff (p, r') i^(q, t") = S^(p + q) <W / 2 (P, r>, t") . 

(34) 



The expression of /2(p, t' , r") can be read from 
eq. (13). Since eq. (33) contains time derivatives 
of Ao.(k, t'), the two point-function of the graviton 
will contain terms such as <9 T '/2(q, t', T")9 T '/ 2 (k — 
q, r', r"), that can be simplified by integration by parts 
in t' or t" and by using the property <9 2 ,/2(q, r', r") = 
-^(r')/ 2 (q,r',r"). 

After several such integrations by parts we can write 
the two-point function of the graviton as 



(Mk, T)h x ,{V, r)) = 6 ^3 M Y ] J dr'dr" J d 3 p.f 2 (p, r', r")/ 2 (k-p, r', r") 



(35) 



x { X i(k, p) [<5"(t') G"(r") + (p 2 + (k - p) 2 ) (g"(t') G(t") + G(r') G"(r")) + 2 (p 2 + (k - p) 2 ) 2 G(r') G(r") 
p |k - p| X2 (k, p) [2 (g"(t') G(r") + G(r') G"(r")) + 4 (p 2 + (k - p) 2 ) G(r') G(r")] } , 



+ 



where we have defined 

Xl (k, p)^^|nf(k)e«(p)e^(k-p)| 2 , 

<7<j' 

X2 (k, p)E^,a' |nf(k)e°(p)e>,(k-p)| 2 , 



& { / ) = 9^1 t G"(r')^5 2 / Gfc(r ' r,) 



a(r') 2 



a(r') 



A2 



(36) 



Up to this point all the expressions used are exact. 
Next, we employ the same approximations of section 
2 A above: we impose w/s(r') ~ m(r') and we drop 
all terms contain fast-oscillating phases of the form 



e 2»/m(f)df_ In this regime 



/ 2 (p, r', r")/ 2 (k-p, r', r") 



(37) 



4m(r') m(r") 

[|/?(p)| 2 |/3(k - p)| 2 + Re [a(p)a*(k - p)/T(p)/3(k - p)]] 



Now all integrals appearing in eq. (35) can be com- 
puted explicitly. In the limit etp/H 2 3> 1 the terms 
in G"(t') turn out to be subdominant with respect to 
those in G(r'). Also, the integral J dr' 'G(r') / 'm(r') is 
suppressed by powers of /ct for fcr 3> 1, while the 
integral in d 3 p is dominated by pr <~ eip/H 2 ^> 1, so 
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that we can set |k| -C |p|. The final result is 



A. The amplitude of the tensor modes. 



dr 



, G(r') 
m(r') 



d 3 pp 4 ( X i(k, p) + X2 (k, p)) (|/3( P )| 2 + 2 |/3(p)| 4 ) . 

(38) 



The integral in dr' appearing in the equation above has 
the same form as that appearing in eq. (18). Therefore 
we treat it as we did in eq. (20). Also, in the regime 
|k| <C |p| we are interested in, xi(k, p) + X2(k, p) = 
g (l — (pk) 2 /p 2 k 2 ) , so that the two-point function of 
the graviton, once we sum on the two helicity states A = 
±1 amounts to twice as much as that of the gravitons 
produced by the scalar sources (19): 



AA' 



k 3 



TT 5 M p 



1 



4^2 



. v 3/2 2 

gip\ (fcr cos kr — sin fcro) 




k 3 



■ ; ■ (39) 



We conclude that the amplitude of the gravitational 
waves induced by vectors equals, per helicity mode, that 
induced by scalars and is therefore negligible. 



The production of tensor modes by a pseudoscalar 
inflaton through gauge field production was discussed 
in [20] and, in greater detail, in [22], where it was 
pointed out that these modes are chiral. We sketch 
here the derivation of the spectrum of gravitational 
waves generated by this mechanism, referring the reader 
to [22] for a thorough discussion. 

In terms of the vector potential A (r, x), defined by 
a 2 B = V x A, a 2 E = —A' and neglecting the spatial 
gradients of <f>, the equations for the gauge field subject 
to the coupling (40) read 



V-A = 0, (41) 



where the prime denotes differentiation with respect to 
the conformal time r. We promote the classical field 
A(t, x) to an operator A (r, x), which we decompose 
on a basis of helicity vectors e± as in eq. (32) . 

The functions A± must satisfy the equations A'± + 
(k 2 =F k<j)'/f)A± = 0. Since we are working on an 
inflating background, we assume the de Sitter met- 
ric a(r) ~ -1/(Ht) and <f> = cp'/a = V2eHM P ~ 
constant. Hence, the equation for A± reads 

d2A± d ^ k) + (k 2 ±2k^A ± ( Tl k) = 0, (42) 

where we have defined 



4. GRAVITATIONAL WAVES PRODUCED 
DURING AXION INFLATION VIA HELICAL 
PHOTONS. 



In this final section we study a scenario where a pseu- 
doscalar inflaton <f) interacts with a gauge field F^ v 
through the coupling 



La 



4/ 



(40) 



where / is a constant with the dimension of mass. The 
rolling inflaton excites, through this coupling, quanta 
of the electromagnetic field, which in their turn source 
the tensor components of the metric. 

In this section we will first derive the amplitude and 
the properties of the spectrum of tensor modes gen- 
erated by this mechanism. We will then study the 
prospects of direct detection of such modes, focusing 
on the specific case where V(<p) oc 4> 2 . 



2fH 



e M P 
2 T 



(43) 



We will be interested in the case £ > 0(1), and we 
assume, without loss of generality, that £ > 0. Then 
A_ stays essentially in vacuum, and we will ignore it 
from now on. However, the mode function A + develops 
an instability, and it peaks at momenta k for which 
(8£) -1 <C \kr\ <§C 2£, where it is well approximated by 



A+(t, k) 



1 



2k \2£aH 



1/4 



(44) 



The "wrong" sign of the term proportional to £ in 
eq. (42) induces an exponential amplification oc e 71 ^ of 
the mode function for the photon A + at sub-horizon 
scales. Such a large occupation number for the vec- 
tor field is in its turn a strong source of gravitational 
waves. Another inflaton-gauge field interaction term 
that leads to amplification of gauge field modes is the 
kinetic coupling f(<p) F^ v . Such coupling, however, 
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generates only a moderate occupation number of su- 
perhorizon photons [29] and, as a consequence, is not 
expected to induce an important production of tensors. 

We can now study the production of gravitational 
waves induced by the helical photons. To do so, we 
plug the expression (44) into eq. (33), and we project 
onto its left- and right-handed components. Of course, 
one should also take into account the parity-symmetric 
component of gravitons that is generated by the usual 
amplification of vacuum fluctuations in de Sitter space. 
This is uncorrelated from those discussed above so that 
the overall power spectra of the hclicity-± components 
of the graviton can be written for £ > 2 (as we will see, 
this is the regime we are interested in) as 



V 



t,+ 



H 



ir 2 M 2 



1 



.6 x 10 



- 7 h 2 



Ml 



H 2 



tt 2 M 2 



1 + 1.8 x 10~ 9 



H 2 e 4 -« 



M 2 



(45) 



We thus see that, as a consequence of the violation of 
parity, the amplitude of the spectra of the left- and 
the right-handed tensor modes generated by the gauge 
field differs by a factor ~ 10 3 . While the parity vi- 
olating component could in principle be exponentially 
large, it was pointed out in [20] that the gauge field 
also contributes, through its coupling with the inflaton, 
to the spectrum of scalar perturbations. This contribu- 
tion is highly nongaussian, and its amplitude is there- 
fore strongly constrained by the non-observation of non- 
gaussianities in the Cosmic Microwave Background and 
the Large Scale Structure bispectra. It turns out that 
the parameter £, when computed at CMB and at LSS 
scales, is constrained to be smaller than about 2.6. This 
implies that the ^-dependent contribution to the tensor 
spectra (45), when computed at cosmological scales, is 
negligibly small. 

Now, the main observation of this section is that the 
quantity £ = 4>/(2 f H) is time dependent and increases 
as the inflaton rolls down its potential. The condition 
£ ;$ 2.6 originates from constraints from CMB and LSS 
data. Therefore £ had to be smaller than 2.6 when 
LSS scales exited the horizon. However, £ can be much 
larger at later times when scales relevant to gravita- 
tional wave interferometers left the horizon. Since there 
are some 40 efoldings of inflation between the time LSS 
scales left the horizon and the time LIGO scales left the 
horizon [8], it is necessary to consider the entire shape 
of the inflationary potential to know how £ evolves. In 
the next subsection we will study in detail, as an exam- 
ple, the parameter space for this scenario in the case of 
a quadratic inflationary potential. As we will see, there 



10" 16 10" 12 10" 8 10 



f(Hz) 



FIG. 1. Amplitude of gravitational waves as a function of 
frequency in the case A4 = 55, £c = 2.1. The star denotes 
the projected sensitivity of Advanced LIGO. 



exists a portion of parameter space where Advanced 
LIGO will be able to observe gravitational waves even 
if the current bound from nongaussianities is satisfied. 



B. Gravitational waves from natural chaotic 
inflation observable by Advanced LIGO. 

In this subsection we study the power spectrum of the 
gravitational waves (45) in the case where the inflaton 
potential takes the chaotic form V(4>) = /i 2 2 /2. The 
model of natural chaotic inflation of [30] (see also [31]) 
leads precisely to this situation: a pseudoscalar inflaton 
with a quadratic potential. Other forms of the poten- 
tial for a pseudoscalar inflaton were considered in the 
models of inflation from axion monodromy [32-34] , and 
we expect that the predictions from these models will 
not differ significantly from those presented here. 



In the case of the chaotic potential V((f>) = fi 2 (j) 2 /2, 
the value of </> during inflation is related to the num- 
ber Af of efoldings before the end of inflation through 
4>{Af) — 2\/U Mp. The parameter /i can be de- 
termined by COBE normalization. Denoting by A/'c 
the number of efoldings corresponding to COBE scales 
(47 < Afc ;$ 62 depending on the details of reheat- 
ing [35, 36]), we have 



6tt 2 P, 



M P , 



(46) 



where P c ~ 2.5 x 10~ 9 . 
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Using the slow roll parameter e = 1/(2 Af), we can 
write £ = y/e/2M P /f as £(A0 = £ c where 
£c is the value of £ computed Afc efoldings before the 
end of inflation. We require £c fi 2.6 in order not to gen- 
erate nongaussianities in the CMB which are too large. 
Note that a comparable constraint on nongaussianities 
originates from the large scale structures at wavenum- 
bers k ~ 0.1 Mpc -1 (see [37] for a recent review), which 
correspond to scales that left the horizon about Afc — 5 
efoldings before the end of inflation. 

Inserting the above expressions into eq. (45), we ob- 
tain the following expression of the energy density in 
gravitational waves as a function of Af 



GW ' 



6 x 10 



-14 



Af c 



(47) 



1 + 4.2 x 10" 



Since gravitational waves at a frequency / correspond 
to scales that exited the horizon about 35+log( //0.1Hz) 
efoldings after the COBE scales, we can plot View h 2 
as a function of / for given Afc and £c by setting 
N = Af c - 35 - log(//0.1Hz) in equation (47). We 
plot in figure 2 the spectrum of gravitational waves for 
a representative set of parameters. 



1. Constraints and detectability. 

The parameter space of this system is constrained by 
the following requirements: first, the backreaction of 
the electromagnetic modes on the inflating background 
must be negligible; second, the nongaussianities induced 
by the same electromagnetic modes through the mecha- 
nism discussed in [20] must be within the limits imposed 
both by CMB and by LSS observations. 

Backreaction on the inflating background is negligible 



for 



/e < 700 V'(<f>) z /H b 



i.e. 



Afc) 



< 6 x 10 



10 



(48) 



This implies that, if we want backreaction to be negli- 
gible all the way to a frequency /, this condition must 
be satisfied with Af = Afc - 35 - log(//0.1Hz). 

As for nongaussianities, the bound of [20] , when eval- 
uated at COBE scales, gives the constraint £c < 2.6. 
Comparable bounds also apply to nongaussianities eval- 
uated at Large Scale Structure scales that left the hori- 
zon some 5 efoldings after COBE scales. We therefore 
impose £(A/"c - 5) < 2.6. 




Nc 



FIG. 2. Values of £c corresponding to detectable tensor 
modes by Advanced LIGO, as a function of the total num- 
ber of efoldings of inflation from the time COBE scales left 
the horizon. The shaded area on the top left corner corre- 
sponds to the region where backreaction cannot be neglected 
and our analysis cannot be trusted. The shaded area on the 
top part of the plot corresponds to the region where LSS 
nongaussianities are too large to be consistent with obser- 
vations. The shaded area at the bottom corresponds to the 
region where the amplitude of tensor modes is below the 
Advanced LIGO detection threshold. Finally, the thinner 
dotted line corresponds to the lower limit of portion of pa- 
rameter space accessible to an instrument such as the Ein- 
stein Telescope. 



We can now discuss the detectability of the tensor 
modes (45) by gravitational interferometers. To fix 
ideas we will focus on Advanced LIGO, which will start 
taking data in the next few years. Advanced LIGO 
is expected to be able to detect flew h 2 = 10~ 9 at a 
frequency of about 100 Hz. The white area in figure 3 
corresponds to the region of parameter space where pri- 
mordial tensor modes might be detected by Advanced 
LIGO without contradicting the constraints described 
above. Detection would be possible for values of £c of 
the order of 2, corresponding to / ~ 10 17 GeV. The 
shaded area in the upper-left corner of the figure corre- 
sponds to a region of the parameter space where back- 
reaction of the electromagnetic modes on the inflating 
background cannot be neglected, and an analysis simi- 
lar to that of [38] is needed. While such an analysis is 
beyond the scope of the present work, it is worth stress- 
ing that that region cannot be excluded by existing data 
and might lead to detectable tensor modes. 

An instrument such as the Einstein Telescope [7] 
would be a factor ~ 10 2 more sensitive in energy than 
Advanced LIGO while working at the same frequencies. 
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The thin dotted line in figure 3 delimits the region of pa- 
rameter space that would lead to a detection of tensors 
by such an instrument. 

Space-based interferometers like LISA, which are sen- 
sitive to much lower frequencies, will not be able to de- 
tect the tensors (45). Indeed, LISA scales are too close 
to cosmological scales and the evolution of £ during in- 
flation is not sufficient to overcome the constraints from 
nongaussianities . 

We also note that the gravitational waves produced 
this way will be chiral [22] and that [39] discussed the 
prospects of a direct detection of a background of chi- 
ral stochastic gravitational waves. These gravitational 
waves will also have a large [11], parity violating [13, 14] 
three point function. Finally, it is nice to speculate how 
a detection of chiral gravitational waves such as those 
described by (45) could correlate with the detection of 
nongaussianities such as those discussed by [20] at cos- 
mological scales: observation of both the nongaussian 
signal in the CMB and of gravitational waves at LIGO 
scales would provide a test of inflation at very different 
times. 

To summarize, we have shown that particle produc- 
tion during inflation generally leads to a feature in the 
spectrum of tensor modes. If such a feature happens 
at the right wavelengths, the tensors might be directly 



detectable by gravitational interferometers. We have 
seen that production of scalar or vector quanta during 
inflation generates only a modest amount gravitational 
waves. However, the analysis of section 4 shows that, if 
the inflaton is an axion coupled to a U(l) gauge field, 
it can lead to an amplification of gravitational waves 
that would be directly detectable by Advanced LIGO 
and/or Advanced Virgo in the next few years. 
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